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Abstract 

Phase transformation in finite-size ferroelectrics is of fundamental relevance for understand- 
ing collective behaviors and balance of competing interactions in low-dimensional systems. We 
report a first-principles effective Hamiltonian study of vortex-to-polarization transformation in 
Pb(Zro.5Tio.5)03 nanoparticles, caused by homogeneous electric fields normal to the vortex plane. 
The transformation is shown to (1) follow an unusual macroscopic path that is symmetry non- 
conforming and characterized by the occurrence of a previously unknown structure as the bridging 
phase; (2) lead to the discovery of a striking collective phenomenon, revealing how ferroelectric 
vortex is annihilated microscopically. Interactions underlying these behaviors are discussed. 

PACS numbers: 77.80.Bh, 77.22.Ej, 64.70.Nd, 77.84.Dy 
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Vortex is a circularly ordered structure phase of broad interest, which has been found 
in ferrous (FM) fi, s „d U c tors g aad Bo S e-E illst ein condensed a toms g 
Recently, the vortex phase formed by electric dipoles, as hinted by Ginzburg et al. when 
exploring the possibility of superdiamagnetism [4], was revealed to exist in ferroelectric (FE) 
nanoparticles. [5| It was shown [5j that the macroscopic toroid moment undergoes a sponta- 
neous transition from being null at high temperature to nonzero at low temperature, acting 
as the order parameter of the vortex phase. In FEs, another ordered structure of general 
concern is the ferroelectric phase with an out- . Transformation be- 

tween these two ordered structures in FE nanoparticles, by electric fields perpendicular to 
the vortex plane, is a phase transition of fundamental importance, since (i) Phase trans- 
formation in finite-size systems has been long known as a fundamental subject [7|; (ii) The 
vortex-to-polarization transformation involves two phases of completely different order pa- 
rameters (namely, toroid moment in vortex phase versus polarization in ferroelectric phase), 
and is a critical example to study collective mechanism in terms of how order parameter 
nucleates and/or disappears within a single domain; (iii) Participation of external electric 
fields during transformation alters and re-organizes the delicate balance among competing 
interactions inside FEs [8], and may lead to new structural phases that are not possible in 
the zero-field case; (iv) It concerns the interaction between FE vortex and electric fields. 
Vortex, as a whole, does not interact with electric fields because of vanished net polariza- 
tion; but individual dipoles do, which explains why collective behavior during the vortex 
transformation is interesting. 

Technologically, the FE vortex phase and the polarization phase both are of considerable 
significance. 9|, [lOj The former promises to increase the storage density of nonvolatile ferro- 
electric random access memories by five orders of magnitude [5], while FEs in polarization 
phase exhibit large strain and electromechanical response, which is of widespread use in 
p—c tra n sd ace re and actual, Q Nowaday FE particlesQ , wires [13fl, and 
nanotubes 14j were all made possible in experiments. 

Despite the fact that vortices exist both in FE and FM particles, they are different in 
key aspects arising from the profound difference between electrostatic and magnetostatic 
interactions. First, the demagnetizing energy in FM particles is 10 3 times smaller than 
the depolarizing counterpart in FE particles of the same size, leading to the fact that for 
particles less than lOnm the dominant interaction in FMs is the short-range (SR) exchange 
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interaction which favors magnetization, much like in bulk 15j, [16| . As size increases, the 
demagnetization energy (scaling approximately as N§ where N is the number of dipoles 
in particle) becomes more important than the SR exchange energy (scaling as N), thereby 
vortex turns favorable until domains start to play a role. Indeed, the minimum threshold size 
for magnetic NigoFeMMos particles to adopt a vortex structure is experimentally found as 
100 nm in diameter (particles less than this size exhibit instead magnetization). [171] This is in 
sharp contrast to the situation in FE nanoparticles where the depolarizing energy is dominant 
and which exhibit a vortex phase as small as 4 nm.[5] The balance between LR and SR 
interactions is thus profoundly different in two materials. Furthermore, because of the large 
depolarizing energy, the shape and anisotropy are insignificant for FE particles when size is 
small (< 10 nm). In fact, FE vortices have been found in particles of cylindrical cubic [3], 
and rectangular shapes. As a result of these differences, the knowledge of transforming 
magnetic vortex is not suitable for FE nanoparticles. 

Here we perform ab-initio based simulations on vortex-to-polarization phase transforma- 
tion, induced by homogeneous electric fields normal to the vortex plane, in nanoparticles 
made of technologically important Pb(Zr . 5^0.5)03 (PZT) solid solution. These studies lead 
to (1) the discovery of a critical new structure that bridges the transformation between two 
phases of unlike order parameters, (2) an unusual collective mechanism showing how FE 
vortex is annihilated, (3) the existence of a previously unknown hysteresis, and (4) an in- 
novative approach for reading FE toroid moment. Neither the phase transformation path 
reported here, nor the unusual dipole behaviors associated with this path, occur for inplane 
electric fields (which cause only dipole flipping [18|). Further, inplane fields do not produce 
out-of-plane polarization that is most utilized in practice. 



We use first-principles derived effective Hamiltonian [19|, [20J and Monte Carlo (MC) 
simulations to determine dipole configuration. Calculations are performed for PZT nanos- 
tructures of cylindrical shape, with diameter d ranging from 9 to 25 and height /i=14 (both 
d and h are in units of bulk pseudo-cubic lattice parameter a=4.0A). The crystallographic 
[001] direction is chosen as the cylindrical 2-axis, with the [100] and [010] directions being the 
x and y axes. Initial dipole configuration of vortex is obtained from annealing simulations 
with temperature decreased in small step. Electric field of varied strength is applied along 
the +z direction at a fixed temperature of 64K, and is coupled with dipole pj at cell i by 
• E ■ pi as described in RefJ2ll. Dipole-dipole interaction in nanoparticles is handled in real 
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space. [3] Typically 10000 MC sweeps are used to simulate dipole responses at each field. 
Outcome of simulation is the soft-mode vector field {uj} (which is directly proportional to 
local dipole p,), homogeneous strain {^}, and inhomogeneous cell shape. 

Figure dh depicts the collective behaviors of toroid moment, G = ^k~Yli r i x P«> anc ^ 
net polarization, P = Y2i Pi (where N c is the number of bulk cells and Q is the cell 
volume), that develop in a d=19 nanodisk as the strength of electric field varies. When 
electric field is small and below a critical value E c l = 1.5V/nm, the disk shows only a G z 
toroid component while G x and are null. The system in the E < E c l field region thus 
retains the same macroscopic toroid symmetry as in zero field (This structure phase will 
be denoted as phase I hereafter). As the field reaches E cl , the G z moment declines only 
slightly as compared to the initial zero-field value, and meanwhile, a net polarization of 
0.3C/m 2 develops. In fact, this mag nitude of polarization is large and comparable to the 



value found in bulk BaTiCh. 
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2l| The vortex response of phase I is thus characterized 



by the coexistence of strong toroid moment and large polarization, both pointing along the 
cylindrical z-axis. We further numerically find that the toroid moment in phase I responds 
to the E field by accurately following a quadratic scaling law as G Z (E) = Go — aE 2 , where 
Go is the zero-field G z moment and the a coefficient is determined to be 285.5 eA 4 /V 2 for 
the d=19 disk. 

The system behaves in a markedly different fashion as the electric field exceeds E c l , man- 
ifested in FigilK by the dramatic occurrence of a nonzero G y component and simultaneously 
a sharp decline of the G z . Being perpendicular to the initial G z moment, the appearance 
of the G y moment deviates the system from continuing to possess macroscopic cylindrical 
symmetry, and the resulting symmetry-broken new structure is to be denoted as phase II. 
Thus for the first time the ferroelectric vortex, under the application of a uniform perpen- 
dicular field, is discovered to be able to self-organize to generate a lateral toroid moment. 
Interestingly, this behavior has never been reported for ferromagnetic vortices (to the best 
of our knowledge). We have also performed calculations for pure PbTiOs nanodisks and 
found that the G^ toroid moment robustly occurs, revealing that the existence of phase II 
is not related to the composition fluctuation in PZT. Assuming that dipoles respond by ro- 
tating collectively toward the field direction, one would anticipate that the drastic decrease 
of the G z moment shall be accompanied by a sharp rise in the polarization. Surprisingly, 
the net P 2 polarization in FigH^, apparently does not feel the drastic variation of the G 2 , 
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and remains, to the naked eye, fairly smooth. As the electric field exceeds a second critical 
value E Cj 2=2.8V/nm, all G components vanish and the system becomes a phase of pure 
polarization. Low-symmetry phase II thus acts as the key intermediate state for bridging 
the transformation from phase I (of the same cylindrical symmetry as the initial vortex) to 
the destination phase of uniform polarization (that also has cylindrical symmetry). 

The evolution of toroid moment during the vortex transformation is also reflected in 
dielectric %33 susceptibility (FigJTb), showing a noticeable hump at E C) i field. Furthermore, 
near the critical E c 2 field where vortex collapses, the X33 coefficient exhibits a sharp decline, 
demonstrating that dipoles in vortex phase generate a notably larger dielectric response than 
dipoles in ferroelectric phase do. 

We now present microscopic understanding of the puzzling collective behaviors that FE 
vortex displays in Figfjji along the transformation path. FigfS] shows snapshots of the 
dipole configurations corresponding to four selected fields Ej (i=l~4) labelled in Figflji. 
Comparison of the dipole patterns at Ei and E2 fields shows that the vortex response in 
phase I is characterized by collective rotation of dipoles towards the field direction. Note 
that at E2 field the vortex pattern on the xy plane maintains a cylindrical symmetry. The 
rotation is thus symmetry conforming in the sense that dipoles within the same distance 
from the central axis respond equivalently. Another important conclusion from simulations 
concerns where the polarization nucleates, since the phase transition involve two different 
order parameters. We find that the z-axis polarization is nucleated at the center of the FE 
vortex. In other words, formation of polarization begins with and is initiated by the rotation 
of those dipoles near the cylindrical axis, which is consistent with the explanation that a 
large strain exists in the vortex center and the rotation reduces this strain energy. 

As the field changes slightly from E2 to E3, two striking phenomena occur: (1) The dipole 
components in the xy plane cease to exhibit cylindrical symmetry (see the E 3 pattern, lower 
panel, Figj2j) . In fact, the inplane components disappear remarkably for dipoles within a 
certain, but not across the entire 360° degree of, azimuthal angle. (2) A lateral G y vortex - 
for which the corresponding toroid moment points at the y-axis — starts to nucleate near 
the right-side surface on the xz plane (see the E 3 pattern, upper panel of Figj2]) . Note that 
this G y vortex forms neither at the cylindrical axis, nor at the same side as the azimuthal 
angle where inplane dipole components are annihilated, in order to avoid the large strain 
energy when vortex develops. 
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Further increase of the electric field from E3 to E4 leads to the widening of the annihilation 
angle (the E4 pattern, lower panel, Figf5]) and full development of the lateral G y vortex. The 
system at the E 4 field is thus characterized by formation of a vortex-free ferroelectric region 
on one side — and a vortex region with perpendicular G y moment on the other side (the E 4 
pattern, upper panel, FigfS]). Interestingly, despite that formation of the lateral G y vortex 
forces some dipoles to point opposite to the field direction between E 3 and E 4 fields, the 
net P z polarization nonetheless remains increasing as seen in FigJTk, consistent with the 
Le Chatelier's principle^. Finally, when electric field continues to increase above E 4 , the 
ferroelectric region expands by extruding the center of the G y vortex toward the right-side 
surface, and this G y vortex eventually disappears at critical E c 2 field. 

Though complex, the dipole behaviors in Figf2]have a simple explanation, namely they 
result from the competition between the applied electric field and depolarizing field. For 
phase II at E c> \ < E < E c 2 , these two fields becomes comparable, and the nanodisk desires 
to reduce the z-axis depolarizing field by forming lateral G y vortex. Indeed, as shown in 
Figj2] for the system at E 4 field, the depolarizing field in phase II is mainly confined to 
the ferroelectric region while being reduced in the vortex region. Meanwhile, to enhance 
the interaction between polarization and the external field, those dipoles in the ferroelectric 
region favor to rotate toward the z-axis more than the dipoles in the G y vortex region, leading 
to the collective phenomenon that the inplane G z vortex is annihilated without maintaining 
the azimuthal symmetry. 

The presence of phase II (with toroid moment rotated by 90° with respect to the initial 
G^ vortex) raises a question of what may happen when one starts with this phase (e.g., at 
E=2V/nm in FigJTk,) and then decreases the external field. Our simulations reveal that, 
regardless of whether the field is gradually reduced or suddenly switched off, the system in 
phase II does not transform back to the initial state of G z moment, and instead is trapped 
in the G y vortex state. More specifically, as the field is reduced, the G y vortex at the E 4 
field in Figj2] grows by moving its core toward the cylindrical axis (not toward the right-side 
surface), resulting in a pure vortex state with a G^ moment and P=0. This leads to a 
hysteresis as shown in the left plot of FigJTb, which is interesting in the sense that (i) the 
hysteresis is caused by toroid moment, not polarization; (ii) it exists in a single particle of 
nanometer size; (iii) during the hysteresis the toroid moment is rotated rather than switched 
to the opposite. We further find that the G y and G z phases at zero field in Fig{lb are very 
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close in energy; the latter is lower by ~1 meV per 5- atom unit cell. Trapping of the system 
in the G y state also suggests that this state is stable and surrounded by energy barrier. To 
confirm this, the system of the G y phase is heated at zero field to a chosen temperature T 
and then cooled down to 64K. We found that only when T is above 500K the G y phase is 
able to overcome the barrier and becomes the G z phase (see the right plot of Fig{Tb). 

Our analysis shows that, in addition to the reduction of the depolarizing field, there 
is another factor that facilitates the transformation from phase I into phase II, that is, the 
interaction between local mode and strain. Fig{TH depicts the strain components at different 
fields. At zero field the lattice of the vortex state is pseudo-tetragonal with cj a ratio less than 
1, i.e., with strain components i] xx = t] yy > i] zz , because all dipoles are lying in the xy plane. 
As the field increases, the c/a rises as a result of the polarization-strain coupling. Notably, 
the field at which c/a becomes 1 (i.e., the system becomes pseudocubic) is very close to the 
critical E c l field where phase I is transformed into phase II. The mode-strain coupling jl9| 
— \B\(rj xx u 2 x + rj yy u 2 y + r] zz u 2 z ) (where B is the coupling strength) advances the transition 
into phase II largely due to the increase in atomic volume, AQ/Q = r\ xx + r\ yy + i] zz , which 
is nearly a constant in phase I and rises sharply for E>E Cj i (FigJTJl). 

We next address how the diameter of nanodisk may influence the vortex transformation. 
Interestingly, we find that there is a critical size d c =17, below which the transformation 
path turns out to be different. Simulation results of a d=9 disk (not shown here) reveal 
that the vortex undergoes a continuous transformation into a single domain ferroelectric 
phase, without the bridging phase II. The dielectric %33 coefficient for this disk at zero field 
is determined to be 55 (Fig{Tb) and is considerably larger than the value of 25 in the d=19 
disk, showing that the same electric field is able to induce a much larger polarization in 
smaller disks. The critical E Cj 2 field in the d=9 disk decreases to 1.3V/nm. One main 
difference between d=9 and d=19 nanodisks is the depolarizing effect which is small in the 
former and thus allows a symmetry conforming transformation. 

Finally we point out that the predicted FE- vortex responses have important technological 
implication. The quadratic law, G Z (E) = Gq — crE 2 , tells us that, when a FE nanodisk is 
exposed to an alternating field E(i) = E cos cut, the toroid moment G(t) will respond with 
a double frequency 2u, and its radiation field can be separated from the field of the vibrating 
polarization which responds only with u. The signal with 2uj frequency may further tell us 
whether it is associated with moment G or — G, since the latter field is phase shifted by 
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22j This may thus open a novel approach by using electromagnetic fields of pulse lasers to 
probe and/or read FE vortex state. Compared to mechanical approach using piezoelectric 
force microscope tips, the optic approach is fast and can be performed simultaneously in a 
large quantity 

In summary, (i) the transformation between vortex and ferroelectric phases in FE particles 
is predicted to follow an unusual and symmetry non-conforming path, which is a macroscopic 
manifestation of the delicate competition of microscopic interactions. A new structure phase, 
with rotated vortex moment and coexisting polarization, was revealed as the intermediate 
phase bridging the transformation, (ii) We discovered a striking collective phenomenon 
leading to the annihilation of FE vortex, namely that the vortex does not disappear shell- 
by-shell, but in a peculiar azimuthal annihilation mode, (iii) The center of FE vortex plays a 
crucial role for the nucleation of polarization, (iv) The existence of an interesting hysteresis, 
caused by toroid moment, may act as a previously unknown channel for energy dissipation 
and dielectric loss, (v) Our simulations further suggest an innovative optic approach for 
reading/probing toroid moment. 
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[1] T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, T. Ono, Science 289, 930 (2000). 
[2] A. A. Abrikosov, Sov. Phys. JETP 9, 1174 (1957). 

[3] K.W. Madison, F. Chevy, W. Wohlleben, J. Dulibard,, Phys. Rev. Lett. 84, 806 (2000). 
[4] V.L. Ginzburg, A. A. Gorbatsevich, Yu.V. Kopayev, and B.A. Volkov, Solid State Commun. 
50, 339 (1984). 

[5] I.I. Naumov, L. Bellaiche, and H. Fu, Nature 432, 737 (2004). 

[6] D.D. Fong, G.B. Stephenson, S.K Streiffer, J. A. Eastman, O. Auciello, P.H. Fuoss, and C. 

Thompson, Science 304, 1650 (2004). 
[7] L.D. Landau and E.M. Lifshitz, Statistical Physics, Part 1, (Pergamon, New York, 1980). 
[8] R.E. Cohen, Nature 358, 136 (1992). 

[9] J.F. Scott, Ferroelectric Memories (Springer, Berlin, 2000). 
[10] M.E. Lines and A.M. Glass, Principles and Applications of Ferroelectrics and Related Materials 



8 



(Clarendon, Oxford, 1979). 
[11] K. Uchino, Pizoelectric Actuators and Ultrasonic Motors (Kluwer Academic, Boston, 1996). 
[12] S. O'Brien, L. Brus, and C.B. Murray, J. Am. Chem. Soc. 123, 12085 (2001). 
[13] W.S. Yun, J.J. Urban, Q. Gu, and H. Park, Nano Lett. 2, 447 (2001). 
[14] F.D. Morrison et al, Rev. Adv. Mater. Sci. 4, 114 (2003). 

[15] F. Jona and G. Shirane, Ferroelecric Crystals (Dover Pub., New York, 1993), p 374. 

[16] E. Yu. Vedmedenko, A. Ghazali, and J.-C. S. Levy, Phys. Rev. B 59, 3329 (1999). 

[17] R.P Cowburn et al, Phys. Rev. Lett. 83, 1042(1999). 

[18] H. Fu and L. Bellaiche, Phys. Rev. Lett. 91, 257601 (2003). 

[19] W. Zhong, D. Vanderbilt, and K. M. Rabe, Phys. Rev. Lett. 73, 1861 (1994). 

[20] L. Bellaiche, A. Garcia, and D. Vanderbilt, Phys. Rev. Lett. 84, 5427 (2000). 

[21] A. Garcia and D. Vanderbilt, Appl. Phys. Lett. 72, 2981 (1998). 

[22] V.M. Dubovik, M.A. Martsenyk, N.M. Martsenyk, J. Magn. Magn. Mater. 145, 211 (1995). 



9 




E (V/nm) T (K) E (V/nm) 



FIG. 1: (a) Toroid moment G (using the left vertical axis) and net polarization P z (using the 
right vertical axis) in a d=19 disk. Arrows on the top horizontal axis indicate the selected electric 
fields for which dipole configurations are analyzed, (b) Dielectric X33 susceptibility in the d=9 disk 
(using the left vertical axis) and in the d=19 disk (using the right vertical axis), (c) Left: hysteresis 
of toroid moment caused by increasing and then decreasing the electric field in d=19 disk; Right: 
transformation of the G y vortex into a G 2 vortex by heating the system to 500K and then cooling 
down, (d) Strain components rj xx , r] yy and rj zz (using the left vertical axis) and volume expansion 
f]xx + Vyy + Vzz (filled symbols, using the right vertical axis) in the d=19 disk. The horizontal axis 
in (a), (b), and (d) is the strength of electric field. 
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FIG. 2: Configurations of local dipoles on the central xy cross section (lower panel) and on the 
central xz cross section (upper panel) in the d=19 disk, at selected fields Ei, E2, E3 and E4 as 
marked by the arrows in FigQJ,. The magnitude of each dipole is enlarged for clarity. 
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